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Abstract 

The  dependence  of  electron  transfer  (ET)  coupling  element.  Jab-  on  the  length  of  rigid-rod-like  systems  consisting  of 
bicyclo[l.l.l]pentane  (BCP),  cubane  (CUB),  and  bicyclo[2.2.2]octane  (BCO)  monomers,  has  been  investigated  with  the 
use  of  ab  initio  Hartree-Fock  (HF)  method  employing  Marcus-Flush  two-state  (TS)  model.  The  value  of  Jab  decreases 
exponentially  with  increase  in  the  number  of  the  cage  units  of  the  ct-bonded  molecules.  The  calculated  decay  constant, 
/ ? ,  shows  good  agreement  with  previously  reported  data.  For  molecular  length  ^  15  A,  the  value  of  JaB  becomes 
negligibly  small,  suggesting  complete  suppression  of  the  through  bond  direct  tunneling  contribution  to  ET  pro¬ 
cess.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 


1.  Introduction 

Long  distance  intramolecular  electron  transfer 
(ET)  constitutes  the  basis  of  many  chemical  and 
biological  phenomena  in  nature.  For  this  reason, 
ET  in  organic  molecules  has  been  the  subject  of 
numerous  theoretical  and  experimental  studies  [1- 
3],  Recently,  there  has  also  been  interest  in  intra¬ 
molecular  ET  due  to  its  potential  application  in 
molecular  electronic/photonic  devices  [3-6],  An 
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important  property  quantifying  the  electronics 
application  of  molecules  is  the  ET  rate,  which  in  a 
two-state  (TS)  theory  [1-3]  is  related  to  the  ET 
coupling  matrix  element,  VAB. 

The  dependence  of  ET  rate  and  the  coupling 
matrix  element,  FAb,  on  the  structural  features, 
such  as  the  geometry,  stereochemistry,  and  effec¬ 
tive  length  of  the  molecule  has  been  the  subject  of 
several  studies  in  the  past  [7-21].  Two  recent  ar¬ 
ticles  by  Newton  [3]  and  Newton  and  Cave  [5] 
present  excellent  reviews  on  the  subject.  In  general, 
electron  transport  in  molecules  follows  radically 
different  mechanisms  from  those  known  in  semi¬ 
conductor  materials  and  is  affected  by  such  factors 
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as  the  symmetry,  geometrical  arrangement,  spatial 
extension,  and  the  separation  between  the  partic¬ 
ipating  one-electron  orbitals.  Therefore,  a  sys¬ 
tematic  study  of  the  factors  affecting 
intramolecular  ET  constitutes  an  essential  first- 
step  toward  developing  the  foundation  of  molec¬ 
ular  electronics. 

In  this  Letter,  we  present  the  results  of  our  ab 
initio  molecular  orbital  (MO)  study  of  the 
length-dependence  of  ET  in  cr-bonded  rigid-rod 
molecules.  Specifically,  the  systems  studied  are 
oligomeric  chains  (CH2-cage„=1 2  3-CH2)  involving 
cr-bonded  bicyclo[l.l.l]pentane  (BCP),  cubane 
(CUB),  and  bicyclo[2.2.2]octane  (BCO)  monomer 
(cage)  units.  These  molecular  rods  are  of  interest  in 
molecular  electronic  circuits  and  photo-induced 
ET  reactions  [22]. 

In  a  recent  study  [23],  the  ET  properties  of  the 
three  monomers,  BCP,  CUB,  and  BCO  were  in¬ 
vestigated  by  the  ab  initio  MO  method  in  the 
premise  of  the  Marcus-Hush  TS  model  [1-3].  The 
three  molecular  systems  exhibited  the  same  quali¬ 
tative  behavior  for  ET  with  respect  to  the  con¬ 
formation  of  the  donor  (D)  and  acceptor  (A) 
moieties,  represented  for  simplicity  by  two  CH2 
groups,  attached  to  the  two  sides  of  the  cage. 
Specifically,  the  value  of  VAB  was  calculated  to  be 
maximum  and  minimum  when  the  planes  of  the 
two  -CH2  groups  were,  respectively,  planar  and 
perpendicular  to  each  other.  Based  on  this  study, 
an  operational  mechanism  for  a  conformational 
molecular  switch  was  proposed  [23].  Quantita¬ 
tively,  the  BCO  system  was  calculated  to  have  the 
highest  ‘molecular  resistance’  (lowest  magnitude  of 
Fab)  among  the  three  systems,  followed  by  CUB 
and  BCP  in  that  order.  However,  it  was  not  clear  if 
the  same  trend  also  holds  in  the  oligomeric  forms 
of  the  three  systems. 

Previous  studies  have  shown  that  the  ET  rate 
between  two  active  centers  generally  falls  off  ex¬ 
ponentially  with  the  increase  in  their  separation 
length  [14].  This  exponential  decrease  of  ET  rate  is 
characterised  by  a  decay  constant  [S,  whose  value 
depends  upon  the  molecular  system.  Some  repre¬ 
sentative  values  for  fi  derived  from  experimental 
measurements  are:  <0.2-1. 4  A-1  for  DNA  [15], 
1.0-1. 4  A-1  for  proteins  [16],  0. 8-1.0  A-1  for  sat¬ 
urated  hydrocarbons  [9,17,18],  and  0.2-0. 6  A-1 


[19,20]  for  unsaturated  organic  systems  such  as 
phenylene,  polyene,  polyyne  and  phenylenevinyl- 
ene.  Recent  experiments  [24]  give  a  somewhat 
narrow  range  (0.8-1. 2  A-1)  for  /)  in  the  case  of 
DNA  molecule.  It  is,  however,  clear  from  the 
widely  varying  range  in  the  value  of  [>  that  the  ET 
decay  rate  strongly  depends  on  the  chemical  na¬ 
ture  and  electronic  structure  of  the  molecules. 

The  length-dependence  of  ET  in  BCP,  CUB, 
and  BCO  rods  has  been  the  subject  of  a  number  of 
previous  investigations.  Liang  and  Newton  [11] 
have  calculated  /(  for  BCP  and  CUB  with  the  use 
of  ab  initio  Hartree-Fock  (HF)  method  and  Ko- 
opman’s  theorem  (KT)  approximation.  A  semi- 
empirical  calculation  of  the  length-dependence  of 
ET  rate  in  systems  involving  a  BCO  spacer  has 
been  reported  recently  by  Do  Monte  [22].  The 
decay  constant,  [1,  and  the  ET  coupling  matrix 
element,  VAB,  in  systems  involving  CUB  and  BCO 
as  the  spacer  have  also  been  reported  from  rate 
constant  measurements  [9,18,21]  in  organic  sol¬ 
vents.  However,  a  number  of  inconsistencies 
among  the  reported  results  for  VAB  exist  in  the 
literature.  For  example,  a  fluorescence  measure¬ 
ment  study  [18]  on  organic  molecules  containing 
BCO  oligomers  as  the  spacer  shows  a  monotonic 
decrease  in  VAB  with  length.  In  contrast,  ab  initio 
HF  calculations  by  Liang  and  Newton  [11]  using 
KT  approach  show  a  non-monotonic  variation  of 
F\b  with  the  increase  in  the  number  of  BCO  cage 
units  in  a  molecular  rod.  Furthermore,  the  value  of 
VAB  for  CUB  monomer  derived  from  experimental 
rate  constant  measurements  [9]  differs  from  that 
obtained  from  ab  initio  HF  calculations  [23]  by  an 
order  of  magnitude. 

The  present  study  is  an  attempt  to  resolve  the 
outstanding  inconsistencies  in  the  previously  re¬ 
ported  data  and  also  to  develop  a  more  detailed 
understanding  of  the  qualitative  behavior  of  the 
length-dependence  of  ET  in  cr-bonded  systems 
from  a  reliable  and  uniform  theoretical  treatment. 
Results  are  presented  on  the  dependence  of  ET  on 
the  number  of  the  cage  units  between  D  and  A 
moieties,  represented,  for  simplicity,  by  -CH2 
groups.  The  effect  of  the  one-electron  basis  set  on 
the  calculated  ET  coupling  is  examined.  The  cal¬ 
culated  results  are  compared  with  previously  re¬ 
ported  theoretical  and  experimental  data. 
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2.  Theory  and  calculations 

For  the  weak  coupling  between  donor  and  ac¬ 
ceptor  centers,  the  ET  rate  constant  is  related  to 
ET  coupling  matrix  element  VAb  as  [1-3] 

2n 

Ke,  =  y|EAB|-FCWD,  (1) 

where  FCWD  represents  the  Franck-Condon- 
weighted  density  of  states,  which  reflects  the  nu¬ 
clear  modes  of  the  system.  The  electronic  coupling 
element,  VAB,  represents  the  strength  of  interaction 
between  localized  states  A  and  B  participating  in 
the  intramolecular  ET.  Since  FCWD  in  Eq.  (1) 
only  weakly  depends  on  the  D  -  A  separation  in  a 
molecule  [10,18,22],  the  distance  dependence  of  ET 
rate  is  primarily  determined  by  the  distance  de¬ 
pendence  of  electronic  coupling  element,  VAB,  gi¬ 
ven  as  [22] 

I  Eab|2  =  K0  exp(— /jy?DA),  (2) 

where  K0  is  a  kinetic  prefactor,  RDA  is  the  separa¬ 
tion  between  two  active  reaction  centers,  and  [!  is 
the  decay  constant,  a  characteristic  of  the  spacer 
system,  which  determines  the  distance  dependence 
of  ET  rate. 

In  the  present  study,  we  used  the  Marcus-Hush 
TS  model  [1-3]  as  formulated  by  Farazdel  and 
Dupuis  [25]  to  calculate  the  electronic  coupling 
element,  VAB.  A  brief  description  of  the  method  has 
been  presented  in  an  earlier  publication  [23].  In  this 
approach,  the  VAB  is  calculated  directly  from  the  ab 
initio  self-consistent  wavefunctions  of  the  two 
participating  electronic  states,  A  and  B.  In  contrast 
to  the  more  commonly  used  and  somewhat  less 
expensive  KT  approach,  the  present  approach  of¬ 
fers  fully  relaxed  charge-localized  electronic  states 
participating  in  the  ET  process  and  thus  includes, 
to  some  extent,  electron  correlation  effect  in  the 
calculation.  Detailed  descriptions  of  and  compari¬ 
son  among  various  theoretical  methods  used  for 
calculating  VAB  has  been  given  by  Newton  [3]  in  a 
recent  review.  The  details  of  the  formulation  and 
computational  algorithm  of  the  ab  initio  TS 
method  used  here  can  be  found  in  the  original  pa¬ 
pers  by  Farazdel  et  al.  [6,25].  Calculations  of  VAB 
were  performed  with  the  use  of  the  HONDO-8 
electronic  structure  code  [26]. 


In  the  KT  approach,  the  value  of  VAB  for  cat¬ 
ionic  systems  can  be  simply  approximated  as  the 
half  of  the  splitting  between  the  highest  occupied 
molecular  orbital  (HOMO)  and  the  HOMO-1  en¬ 
ergy  levels  in  the  triplet  state  of  the  neutral  mole¬ 
cule.  Due  to  this  simplicity,  the  KT  approach  has 
been  widely  used  for  estimating  VAB  in  complex 
organic  molecular  systems  [3,5,7,10-13].  In  the 
present  work,  we  have  also  estimated  VAB  using  the 
KT  approach  to  compare  with  our  results  ob¬ 
tained  from  a  more  rigorous  self-consistent  meth¬ 
od  [6,25], 

Three  basis  sets,  namely  STO-3G,  3-2 1G  and 
double  zeta-plus  polarization  (DZP)  Cartesian 
Gaussian  functions  were  used  in  the  calculation  on 
the  monomer  units.  For  the  oligomeric  systems 
(n  >  1),  we  could  not  use  the  polarized  basis  set 
(DZP)  due  to  computational  limitation. 

Initially,  the  geometry  of  neutral  rod  structures 
containing  1,  2,  and  3  monomer  units  terminated 
by  -CH3  groups  at  the  two  ends  were  optimized 
with  the  use  of  STO-3G  basis  set.  Subsequently, 
the  -CH3  groups  at  the  two  ends  were  replaced 
by  -CH2  groups  to  model  the  molecular  inter¬ 
connect.  In  all  the  oligomeric  systems  (n  >  1),  a 
staggered  conformation  of  the  a-bonded  cage 
units  was  found  to  be  the  most  stable  structure. 
The  largest  oligomer,  with  n  =  3,  for  each  system 
optimized  in  this  work  is  shown  in  Fig.  1.  The 
optimized  inter-cage  distances  for  the  rods  con¬ 
taining  BCP,  CUB,  and  BCO  are  1.514,  1.511, 
and  1.590  A,  respectively.  The  optimized  geo¬ 
metrical  parameters  for  the  monomer  units  have 
been  reported  in  an  earlier  publication  [23]  and 
did  not  change  in  the  present  calculations.  The 
geometry  optimization  calculations  were  per¬ 
formed  with  the  use  of  the  Gaussian  94  [27] 
electronic  structure  code. 

In  order  to  calculate  VAB  an  electron  was  re¬ 
moved  from  the  neutral  system  to  create  an  open- 
shell  system  with  one  unpaired  electron.  The 
charge-localized  structures  and  wave  functions 
were  obtained  by  varying  the  C(CH2)-C(cage) 
bond  distances  while  keeping  the  cage  structure 
fixed.  The  charge-localized  wavefunctions  were 
used  as  the  initial  guess  for  the  calculation  of  ET 
matrix  element  at  an  average  geometry  of  the  two 
localized  structures. 
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Fig.  1.  (a)  Optimized  CH2-(C5H6)n=,-CH2  structure;  (b)  op¬ 
timized  CH2-(C8Hf,)„=3-CH2  structure;  (c)  optimized 
CH2-(CgFI12)„=3-CFl2  structure.  The  large  (dark)  spheres  rep¬ 
resent  C  atoms  and  the  small  (light)  spheres  are  H  atoms. 

3.  Results  and  discussion 

The  calculated  results  for  VAB  are  presented  in 
Table  1,  together  with  the  available  theoretical  and 
experimental  results  reported  earlier.  The  two  ba¬ 
sis  sets  (STO-3G  and  3-2 1G)  give  very  similar  re¬ 
sults  although  the  magnitude  of  the  VAB  is  larger 
for  the  latter.  For  the  CUB  system,  improvement 
of  the  basis  set  appears  to  have  no  noticeable  effect 
on  Fab-  The  KT  method  also  shows  a  similar  effect 
with  respect  to  basis  set.  The  only  exception  is 
BCO  system,  in  which  the  TS  self-consistent  ap¬ 
proach  shows  an  increase,  while  the  KT  method 
shows  a  decrease  in  the  value  of  VAB  with  the  im¬ 
provement  of  the  basis  set.  For  the  BCP  and  CUB 
rods,  the  results  obtained  for  VAB  from  the  TS  and 
the  KT  approaches  reasonably  agree  with  each 
other.  In  the  case  of  the  BCO  rod,  the  TS  ap¬ 
proach  shows  a  monotonic  decrease  in  the  ET 
coupling  with  increase  in  the  number  ( n )  of  the 
cage  units.  In  contrast,  the  KT  approach  shows  an 
increase  in  the  value  of  Iab  in  going  from  n  =  1  to 
2  and  then  a  decrease  from  n  =  2  to  3  in  the  BCO 
rod.  A  similar  trend  in  the  change  in  ET  coupling 


with  increase  in  the  number  of  the  BCO  cage  unit 
was  also  noted  in  the  KT  calculations  by  Liang 
and  Newton  [11].  Experimental  measurements  [18] 
on  the  rate  constant,  however,  suggested  a 
monotonic  decrease  in  ET  coupling  with  increase 
in  rod  length  of  the  BCO  system.  A  closer  exam¬ 
ination  of  the  results  in  Table  1  suggests  that  for 
molecular  rod  with  n  =  2  and  3,  the  KT  approach 
tends  to  overestimate  the  value  of  VAB.  To  the  best 
of  our  knowledge,  no  experimental  data  on  the 
length-dependence  of  VAB  in  BCP  rods  is  available 
in  the  literature  for  comparison. 

One  notes  (Table  1)  a  difference  of  an  order  of 
magnitude  between  the  value  of  VAB  for  the  CUB 
monomer  calculated  in  the  present  study  by  the  TS 
self-consistent  method  and  the  experimental  value 
estimated  from  rate  constant  measurements  in 
solution  [9].  The  TS  results  show  a  decrease  in  the 
value  of  Fab  by  an  order  magnitude  as  the  number 
of  the  cage  unit  increases  from  1  to  2.  The  corre¬ 
sponding  change  in  the  measured  value  is  much 
smaller  by  only  a  factor  of  2.  However,  in  a  CUB 
rod,  for  n  =  2  and  3,  the  calculated  TS  results  for 
VAB  show  excellent  agreement  with  the  corre¬ 
sponding  experimental  values  estimated  from  rate 
constant  measurements.  This  could  very  well  be 
accidental  since  the  calculations  have  been  per¬ 
formed  on  isolated  systems  while  experiments  were 
performed  in  solution. 

The  optimized  distance  between  the  two  C  at¬ 
oms  of  the  -CH2  end  groups,  where  the  charge 
localization  occurs,  was  taken  as  the  effective 
tunneling  length,  L,  for  ET.  The  calculated  value 
of  L  in  the  systems  investigated  in  the  present 
study  is  listed  in  the  second  column  of  Table  1.  The 
decay  constant  [i  was  estimated  from  an  expo¬ 
nential  fit  of  Fab  with  L.  A  plot  of  ln(|  Fab | )  as  a 
function  of  L  for  BCP,  CUB,  and  BCO  rod  is 
shown  in  Figs.  2-4,  respectively.  The  plotted 
curves  for  BCP  (Fig.  2)  and  CUB  (Fig.  3)  exhibit 
an  exponential  decrease  in  the  value  of  VAB  with 
increase  in  L.  The  same  is  not  true  for  BCO,  where 
the  plots  deviate  from  a  straight  line,  suggesting  a 
somewhat  non-exponential,  but  monotonic  de¬ 
crease  (Table  1)  in  decay  rate  with  increase  in  the 
value  of  L. 

The  values  of  /(  estimated  from  the  calculated 
slope  of  the  best-fit  straight  line  for  the  three  sys- 
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Table  1 

ET  coupling,  \VAB\  (kj/mol),  as  a  function  of  the  distance  ( L )  between  the  donor  and  acceptor  groups  in  molecular  rods  (CET-[cage],,=[_3-CE[+)  involving  BCP,  CUB  and 
BCO  monomer  units 


(Cage),, 

L* 

(A) 

Present  work  (TS) 

Present  work  (KT) 

Previous  results 

STO-3G 

3-2 1G 

DZP 

STO-3G 

3-2 1G 

DZP 

Theory 

Experimental 

(BCP), 

4.84 

36.16 

47.98 

52.00 

44.38 

58.77 

58.85 

54.08b 

(BCP)2 

8.26 

5.04 

9.64 

11.34 

16.63 

1 7. 12b 

(BCP)3 

11.67 

0.65 

1.77 

2.80 

4.83 

5.62b 

(CUB), 

5.64 

25.65 

25.50 

25.25 

36.66 

39.48 

34.73 

2.51° 

(CUB)2 

9.88 

1.49 

0.80 

9.36 

7.68 

0.49c 

(CUB)3 

14.11 

0.07 

0.06 

2.49 

1.64 

0.05c 

(BCO), 

5.62 

2.29 

4.78 

6.83 

11.58 

2.69 

1.93 

(3.15  ±  3.93)b, 

1.6T 

0.79d,  2. 12e 

(BCO), 

9.96 

0.97 

1.23 

5.80 

5.73 

5.54b 

(BCO); 

14.26 

0.08 

0.14 

1.23 

1.29 

1.63b 

a  Distance  between  the  C  atoms  of  the  two  end  -CET  groups. 
bRef.  [11],  obtained  using  KT  approach  and  super  exchange  analysis. 
c  Ref.  [9],  estimated  from  PR  rate  constant  measurements  in  solution. 

JRef.  [13],  obtained  using  KT  approach. 

eRef.  [7],  The  theoretical  results  were  obtained  using  KT  approach  and  the  experimental  results  were  obtained  from  pulse  radiolysis  (PR)  rate  constant  measurements 
in  solution. 
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Fig.  2.  Electronic  coupling  matrix  element,  (kj/mol),  in 
natural  logarithmic  scale,  as  a  function  of  the  effective  length 
(A)  of  the  BCP-rod. 


Fig.  3.  Electronic  coupling  matrix  element,  Tab  (kj/mol),  in 
natural  logarithmic  scale,  as  a  function  of  the  effective  length 
(A)  of  the  CUB-rod. 


terns  are  listed  in  Table  2.  Also  listed  in  the  table 
are  the  available  literature  values  of  /?  for  the  three 
systems.  The  two  basis  sets,  namely  STO-3G  and 
3-2 1G,  used  for  the  calculation  on  the  oligomeric 
systems  (n  >  1),  give  very  similar  values  for  f3.  For 
the  CUB  system,  the  calculated  [3  is  found  to  be 
somewhat  insensitive  to  the  basis  set.  A  compari¬ 


L(A) 


Fig.  4.  Electronic  coupling  matrix  element,  (kJ/mol),  in 
natural  logarithmic  scale,  as  a  function  of  the  effective  length 
(A)  of  the  BCO-rod. 


son  of  the  present  TS  self-consistent  results  of  (3  for 
BCP  and  BCO  with  the  corresponding  value  ob¬ 
tained  with  the  use  of  the  KT  approach  [11]  sug¬ 
gests  that  the  latter  approach  gives  a  somewhat 
smaller  decay  rate  than  the  former.  Our  calculated 
values  of  (3  are,  however,  in  good  agreement  with 
the  available  experimental  results,  with  the  latter 
being  somewhat  smaller  in  magnitude.  A  small 
quantitative  difference  between  the  value  of  (3  for 
CUB  obtained  in  the  present  study  (TS  self-con¬ 
sistent  approach)  and  that  derived  from  the  ET 
rate  measurements  [9]  reflects  the  difference  in  the 
chemical  and  electronic  environments  around 
the  spacer  systems  in  the  two  studies.  For  example, 
the  experiments  were  performed  in  a  solution  in 
which  the  intermolecular  interaction  between  the 
solvent  and  the  solute  molecules  (the  species  under 
investigation)  may  influence  the  measured  ET  rate. 
In  contrast,  the  present  calculations  have  been 
performed  on  an  isolated  molecule.  Furthermore, 
the  experiments  were  performed  on  a  CUB  system 
with  biphenyl  and  naphthyl  groups  attached  at  the 
two  ends,  whereas  in  our  calculations,  the  rigid- 
rod  systems  were  terminated  by  -CEE  end  groups. 
It  is  also  worth  noting  that  a  fair  degree  of  un¬ 
certainty  exists  in  the  experimental  value  of  FAb  for 
the  monomer  estimated  from  rate  constant  mea- 
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Table  2 

» -l 

Calculated  and  measured  values  of  the  ET  decay  constant,  fl  (A  ),  for  BCP,  CUB,  and  BCO  rods 


System 

Present  work 

Previous  results 

STO-3G 

3-2 1G 

Theory3 

Experimental 

BCP-rod 

1.18 

0.97 

0.66 

CUB-rod 

1.40 

1.42 

0.90b 

BCO-rod 

0.79 

0.82 

0.57 

>  0.70° 

“Ref.  [11],  obtained  from  VAB  values  estimated  using  KT  method. 

b  Ref.  [9],  obtained  from  pulse  radiolysis  measurements. 

cRef.  [18],  obtained  from  picosecond  fluorescence  measurements. 


surement  [9].  In  view  of  the  noted  differences  in  the 
chemical  and  electronic  environments  around  the 
CUB  rod  and  uncertainties  in  the  experimental 
result,  the  proximity  between  our  calculated  value 
for  p  and  that  derived  from  experiment  is  quite 
encouraging  and  lends  confidence  to  the  calculated 
results  obtained  in  this  study.  The  calculated  val¬ 
ues  of  p  for  the  three  systems  suggest  that  the  ET 
decay  rate  is  much  stronger  in  CUB  rod  than  those 
in  the  BCP  and  BCO  rods. 

It  is  useful  to  establish  a  quantitative  relation¬ 
ship  between  ET  and  the  effective  length  of  the 
molecular  rods.  An  examination  of  the  results  of 
Table  1  reveals  that  the  ET  matrix  element  in  all 
three  cases  nearly  vanishes  as  the  effective  length 
increases  to  ~14  A.  In  other  words,  in  the  three  a- 
bonded  systems  studied  here,  there  appears  to  be 
practically  no  ET  beyond  a  length  of  ~14  A.  This 
observation  has  important  implications  for  the  ET 
mechanism  in  non-conjugated  systems.  Since  the 
present  theory  of  ET  is  based  on  the  direct  overlap 
of  many-electron  wavefunctions,  it  provides  in¬ 
formation  on  direct  ‘through-bond’  electron 
transport.  Our  calculated  results  of  VAB  suggest 
that  in  the  a-bonded  systems,  the  direct  through- 
bond  tunneling  of  electrons  is  totally  suppressed  at 
a  length  >  15  A.  It  is  interesting  to  note  that  this 
number  is  remarkably  close  to  a  value  of  —20  A 
for  the  maximum  effective  distance  for  direct  tun¬ 
neling  current  measured  for  self-assembled 
monolayers  (SAM)  of  alkane  chains  [28].  This 
property  of  the  rigid  molecular  rods  can  be  use¬ 
fully  exploited  for  the  design  and  fabrication  of 
molecular  capacitors  (MC)  and  molecular  field- 
effect  transistors  (MFET).  Rigid-rod  anchors 
consisting  of  BCP,  CUB,  and  BCO  with 


lengths  ^  15  A  can  be  used  as  ‘insulators’  to 
completely  suppress  direct  tunneling  current  be¬ 
tween  two  active  reaction  centers  in  a  molecular 
circuit. 


4.  Summary  and  conclusion 

We  have  performed  ab  initio  calculations  on  the 
length-dependence  of  the  ET  coupling  matrix  ele¬ 
ment,  T a ti -  in  a-bonded  rigid  molecular  rod 
structures  consisting  of  BCP,  CUB,  and  BCO 
monomer  units.  We  have  estimated  the  ET  decay 
constant,  p,  from  the  calculated  length-dependence 
of  FAb-  Calculations  on  systems  containing  a  single 
cage  unit  were  performed  with  three  one-electron 
basis  sets,  namely  STO-3G,  3-2 1G,  and  DZP.  For 
oligomeric  systems,  the  calculations  were  limited 
to  the  two  small  basis  sets  due  to  the  restrictions 
imposed  by  the  computational  requirements.  For 
the  systems  consisting  of  BCP  and  BCO,  the  cal¬ 
culated  value  of  Fab  exhibits  strong  sensitivity  to 
the  basis  set.  In  contrast,  the  value  of  VAB  in  the 
case  of  CUB  system  appears  to  be  less  sensitive  to 
the  quality  of  the  basis  set.  However,  both  STO- 
3G  and  3-2 1G  basis  sets  give  a  similar  trend  for  the 
distance  dependence  of  VAB,  lending  confidence  to 
our  calculated  values  of  p. 

A  number  of  useful  findings  concerning  ET  in 
a-bonded  rigid  molecular  rods  have  emerged  from 
the  present  study.  Our  calculations  reveal  that  the 
ability  of  a  a-bonded  rigid  molecular  rod  to  act  as 
an  insulator  in  molecular  circuits  depends  upon  its 
electronic  and  geometrical  features.  For  the 
monomeric  units,  the  magnitude  of  F\r  is  calcu¬ 
lated  to  have  the  order,  BCP  >  CUB  >  BCO, 
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suggesting  the  BCO  to  be  the  best  candidate 
among  the  three  for  application  as  a  molecular 
insulator.  However,  with  the  increasing  number  of 
the  monomer  units  in  the  rod,  the  order  in  the 
calculated  magnitude  of  FAb  changes  to 
BCP  >  BCO  >  CUB. 

The  calculated  values  of  VAB  decrease  expo¬ 
nentially  with  an  increase  in  the  effective  length  of 
the  molecule,  which  is  consistent  with  previously 
reported  theoretical  and  experimental  results  in 
organic  systems.  The  value  of  the  ET  rate  decay 
constant,  [!,  estimated  from  the  length-dependence 
of  Fab  has  the  largest  magnitude  in  the  case  of 
CUB  followed  by  BCP  and  BCO  in  that  order. 
This  suggests  that  among  the  three  molecular  rods 
investigated  in  the  present  study,  the  one  consist¬ 
ing  of  CUB  cages  is  more  efficient  in  controlling 
the  ET  rate  and  hence  the  tunneling  current  than 
those  consisting  of  BCP  and  BCO  monomer  units. 
The  />  values  obtained  from  our  calculations 
(0.8-1. 4  A-  )  is  consistent  with  the  experimental 
observation  (0. 8-1.0  A-  )  in  o-bonded  systems. 
Furthermore,  the  present  study  also  lends  support 
to  the  experimental  observation  relating  the  max¬ 
imum  effective  length  of  SAM  of  alkane  chains 
exhibiting  electronic  conduction.  The  present  cal¬ 
culations  reveal  that  the  direct  through-bond  ET 
in  molecular  rods  consisting  of  BCP,  CUB,  and 
BCO  units  is  completely^ suppressed  for  an  effective 
molecular  length  ^  1 5  A.  This  value  can  be  used 
as  a  guide  in  designing  MFETs,  where  the  ct- 
bonded  rods  in  the  molecular  circuit  can  be  used  to 
suppress  leakage  current  from  the  gate  end. 
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